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The CALICE [1^ Si-W electromagnetic calorimeter [2 has been tested with electron 
beams (1 to 6 GeV) at DESY in May 2006, as well as electrons (6 to 45 GeV) and 
hadrons (6 to 80 GeV) at CERN in August and October 2006. Several millions of 
events have been taken at different incident angles (from 0° to 45°) and three beam 
impact positions. The ECAL calibration is performed with muon beams and shows a 
good uniformity for nearly all channels. The large statistics available allows not only 
to characterise the ECAL physics performance, but also to identify subtle hardware 
effects. 



1 Introduction: the Calice ECAL Prototype 

The Si-W ECAL physics prototype is composed of 30 layers of 3 x 3 wafers, each wafer 
having an array of 6 x 6 pixels of 1 x 1 cm^. The two top rows of wafers are completed for 
the full depth in July 2006. The mechanical structure consists of tungsten sheets wrapped 
in carbon fibre, providing 15 alveola where slabs are inserted. One slab is made of two PCBs 
on each side of a tungsten layer, with the wafers conductively glued to the PCB. The very 
front end electronics (VFE) provide preamplification and are located outside the active area, 
but mounted on the same PCB as the silicon wafers. The prototype is built of three stacks, 
each of ten layers of alternating tungsten and silicon. Each stack has a different tungsten 
thickness: 1.4 mm or OAXq per layer in the first stack, 2.8 mm or O.SXq per layer in the 
second stack and 4.2 mm or 1.2Xo per layer in the rear stack. This choice should ensure 
a good resolution at low energy, due to the thin tung- 
sten in the first stack, combined with a good contain- 
ment of the electromagnetic showers, with an overall 
thickness of about 20 cm or 24Xo at normal inci- 
dence. To rotate to angles of 10°, 20°, 30° and 45° 
with respect to the beam axis, the three stacks, me- 
chanically separate, are also translated laterally so 
that the beam still passes through all of them. 

The purpose of the testbeam phase is to validate 
the simulation against a realistic detector, as well 
as allowing to detect the potential hardware prob- 
lems. Once the simulation is trusted, full detector 
studies will lead to the optimisation of the calorime- 
ter for a Particle Flow approach. The test setup at 
CERN is presented on Figure [U and is simulated us- 
ing Mokka [3]. 

Three drift chambers are used for the tracking at CERN and four at DESY. The ECAL, 
HCAL [4j and TCMT [5J follow, with the ECAL mounted on a movable stage for angle 
scans. 




Figure 1: Testbeam setup at CERN, 
in August 2006. The ECAL is 
in front, followed by the analogue 
HCAL and TCMT. 
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2 Summary of the data taken 



Eight million triggers were taken at DESY during 14 days in May 2006, for seven beam 
energies (1, 1.5, 2, 3, 4, 5, 6 GeV), five angles (0°, 10°, 20°, 30°, 45°) and three positions of 
the beam on the ECAL front face (center, border and corner of a wafer), with a minimum 
of 200,000 events per configuration. Six layers were not instrumented: the last eight layers 
had one dummy slab, one instrumented slab, and two dummy slabs. 

The August 2006, CERN beams allowed to take another 8.6 million triggers in ECAL 
only mode, with the full depth instrumented, for six beam energies between 6 and 45 GeV, 
and four angles. Pions were taken between 30 and 80 GeV in combination with the HCAL 
and TCMT. For calibration purposes, 30 million muon events were also taken parasitically 
to an experiment upstream. 

The setup in October 2006 was slightly different, with ECAL and HCAL at 6 cm from 
each other. 3.8 million triggers were taken with electrons and positrons from 6 to 45 GeV, 
and 22 million with pions from 6 to 80 GeV. Another 40 million muon events were added 
for calibration purposes. 



3 ECAL Calibration 
3.1 Gain Calibration 

The current calibration of the ECAL prototype is performed by using a set of 74 high- 
statistics muon runs ( ~ 250,000 events each), taken during October 2006 with another 
experiment upstream, providing a wide spread of the muon beam over the front face of the 
prototype. The events are triggered with a 1 m^ scintillator counter. 



Calibration Constants 




CalibConsts 
Entries 6471 
Mean 45.54 



After pedestal substraction (see Section [3T2|) . 
the runs are reconstructed using a fixed global 
noise cut of half a Minimum Ionising Particle 
(MIP), 1 MIP being estimated to 50 ADC counts 
by former studies. To reject any remaining noise 
hits, it is required that the hits of one event form 
tracks characteristic of a MIP. 

The distributions obtained by channel are de- 
scribed by a convolution of a Gaussian and a 
Landau function. The calibration constant is de- 
fined as the most probable value (MPV) of the 
Landau function, while the standard deviation of 
the Gaussian defines the noise value for each cell. 
For 6403 out of the 6480 channels of the proto- 
type, a calibration constant can be obtained via 
this method without further investigations. The 

remaining cells show a noise value that is unusually high, and are treated by estimating 
the additional noise contribution, or by applying the calibration constant from one of their 
neighbours. 

An entire wafer (36 channels) seems to not be fully depleted at the applied voltage of 
200V, resulting in a MIP peak at half the normal value. A relative value between the mean 



Figure 2: Equivalent ADC value for the 
energy deposit of 1 MIP for the live cells 
of the ECAL prototype. 
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MIP signals of the wafer and its neighbours can be estimated, ahowing a relative calibration 
of the cells. 0.14 % of the channels give no output and were considered as dead. The cali- 
bration constants for all calibrated channels are histogrammed in Figure [21 The distribution 
is narrow, with almost all pads in the range 40 to 50 ADC counts per MIP. The small peak 
at 23.5 corresponds to the single incompletely depleted wafer. 

3.2 Pedestal 

For all beam tests performed, the data acquisition consists of a fixed sequence of 500 pedestal 
events, 500 events with charge injection via the calibration chips, and then 20,000 beam 
events. The pedestal events are used to make a first estimate of the pedestal (mean value) 
and noise (RMS) per channel. 

It has been observed that the pedestals are not necessarily stable, but subject to a 
random shift affecting all channels of one layer with the same drift. This effect concerns 
several particular PCBs, is time dependant, and is attributed to the instabilities of the power 
supplies giving the pedestal lines, which are not isolated. To correct for these instabilities, 
the pedestals are recalculated on an event by event basis, by discarding all cells recording a 
signal, and iterating until the RMS of the distribution obtained with the remaining channels 
is of the order of the expected noise. 

3.3 Noise 



ROOT correlation factor between 2 channels ROOT correlation factor between 2 channels Noise per cell for layer 8 




Figure 3: Layer 8, left: module recording a signal, middle: neighbouring wafer, right: mean 
noise per wafer. Upper row: without corrections, and bottom row: with corrections. 

In order to identify coherent noise, the correlation between pairs of channels is calcu- 
lated in signal events. No clear correlation is seen in an entire PCB, except the one coming 
from the pedestal drift discussed above. The results are thus shown only on a wafer basis. 
Figure [3] displays the correlation factors (colour scale) as a function of the channel indices, 
numbered from to 35, for a particular wafer affected by the pedestal drifts described above, 
and a run recorded at DESY with an energy of 6 GeV. Channels numbered to 17 and 18 to 
35 correspond to two different chips. The corresponding noise level per layer, for all wafers. 
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is also presented in Figure [3l on the right column. 

The results before any corrections are shown in the upper row of Figure [3l The left 
plot represents the module in which the beam was directed. It can be seen that the re- 
gion with signal shows less correlations, due to the fact that most pixels are discarded in 
the noise calculation. The middle plots show a neighbouring module, affected only by the 
global pedestal drift. The difference between these two wafers allows the identification of 
a crosstalk issue. All the channels of some wafers recording a high signal, like the one pre- 
sented in Figure [3] on the left, suffer from a pedestal shift towards negative values. This 
effect does not propagate to the neighbouring wafers. This seems to be random in space 
and in time, but it is clearly correlated with the intensity of the signal recorded. This effect 
is not yet understood, but under investigation. In order to correct for the induced corre- 
lated noise, the mean and standard deviation are calculated on an event-by-event basis, per 
wafer, after discarding the signal hits, and iterating over the channels taken into account in 
the sum. The bottow row on Figure [3] shows the results after all the pedestal corrections 
described above. The corrections are performing well, bringing the noise back to the normal 
level of 6 ADC counts. For this particular layer, affected by both problems (i.e. coherent 
noise on the PCB level due to an instability in the power supply, and crosstalk affecting the 
wafer recording a high signal), the wafer recording a signal still shows a remaining 20% of 
correlations. The correlation is however completely removed for wafers affected only by the 
crosstalk problem, which confirms that the corrections are performing well. 



4 Background for physics performance studies: electron selection 



The selection of the 
single electron events 
is loose, in order to 
avoid bias: 

• Ecell > 0.6 MIP 

removes the noise. 
The threshold 
is about five 
times the aver- 
age noise mea- 
sured per cell. 
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Figure 4: Total ECAL energies for a 3 GeV e" DESY beam {left) 
and for 20 GeV at CERN (right) ^ with the energy selection windows. 
The shaded areas show the effect of the cuts on the shower barycentre 
(left) and on the Cerenkov counter signal (right). 



the total en- 
ergy recorded 
in the ECAL, 
^raw 5 should be 

in the range 125 < (£;raw(MIP)) / (^bea 
three stacks weighted in the ratios 1:2:3 



n(GeV)) < 375. £^raw is computed with the 
according to the tungsten thickness. 



Further cuts are applied to some particular samples: the significant pion content of some 
high energy electron runs at CERN is reduced by using the threshold Cerenkov counter, 
whereas the low energy halo coming with some low energy DESY beams is rejected with 
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additional cuts on the shower bary centre. The effect of these two additional cuts is indicated 
by the shaded regions in Figure [H 



5 Energy Response, Linearity and Resolution 




^15 15000 f i 



dlO 10000 i 




I GAUGE Preliminarj - 



□ m 



The total response of ECAL is computed by summing the hit energies in the three sections 
of the detector. If Ei^ £^2, Es are the recorded energies in the first, second and third stack 
respectively, the total response is £^tot = {o^iEi -\-a2E2 -\- a^E^)/ /3. The naive choice for the 
weights (ai, 0^2, 0^3) = (1,2,3) is generally used. It reflects the relative thicknesses of the 
tungsten layers in each of the stacks, and hence the relative sampling fractions. However, a 
weighting scheme optimisation for energy resolution was performed as well, leading to the 
slightly different values of (1.1, 2., 2.7). The normalisation (3 has been arbitrarily fixed to 
250 MIP/GeV. 

The guard rings create 2 mm non-active inter-wafer gaps, causing non-uniformities in the 
ECAL response, as illustrated in Figure[5l left, where the mean energy is plotted as a function 
of the shower barycentre, (x, y) = ^^{EiXi^ Eiyi) / ^ • Ei. Dips in response are clearly visible 
and account for the asym- 
metric tail on the low 
side of the distribution 
of total energy (Figure O 
right), which is reason- 
ably well modelled by 
the simulation. The 
correction of these non- 
uniformities will be dis- 
cussed in Section 6. 

The beam profile, and 
thus the fraction of beam 
particles traversing the 
inter- wafer gaps, depends 
on the beam energy. 
Therefore, in order to 
avoid bias, a cut is ap- 
plied on the shower barycentre position such as to select showers not affected by the gaps. 

To estimate the en- 
ergy resolution and the 
mean calorimeter response, 
the distribution from Fig- 
ure [6l left, is fitted by 
a Gaussian in the asym- 
metric range of [—a, +2<j] 
in order to reduce sen- 
sitivity to pion back- 
ground and to radiative 
effects upstream of the 
calorimeter, as well as 
to any residual influence 
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Figure 5: Left: mean ECAL energy as a function of the shower 
barycentre (x^y). Right: total ECAL energy for a 30 GeV e~ 
beam, for data (points) and Monte Carlo ( open histogram ). 
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Figure 6: Left: a Gaussian fit to the measured energy, for 30 GeV 
e~ data. Right: ECAL energy response, divided by beam energy, 
as a function of beam energy. 
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of the inter- wafer gaps. 

The ratio of the reconstructed energy to the 
beam energy, as function of the beam energy is 
shown in Figure [6l right, for the two choices of 
weights. Non-hnearities are at the % level. The 
linearity is somewhat better for the optimised 
weights. 

The energy resolution, similar for both 
weightings, is shown in Figure H The 
Monte Carlo prediction, with the naive weights, 
also shown, is in reasonably good agree- 
ment with the data. The resolution can 
be parametrised, for the naive choice of 
weights and, respectively, the optimised one, 
as 





CALICE Preliminary 
• a=(1,2,3) 
o a=(1.1, 2, 2.7) 
^ Monte Carlo a=(1,2,3) 
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Figure 7: ECAL energy resolution as 
function of 1/VE. Two possible weight- 
ings of the ECAL stacks are compared. 
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6 Interwafer gap corrections 
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Figure 8: 15 GeV e", 
rections (blue circles): 




energy (GeV) 



The method used for 
correcting the interwafer 
gaps operates at the 
event level and relies only 
on the calorimeter infor- 
mation: it parametrises 
the mean calorimeter re- 
sponse as function of 
the shower position in 
the calorimeter and ap- 
plies subsequently cor- 
rective factors for each 
event according to this 
parametrisation. It is 
only geometrical, inde- 
pendent in X and y and 

without any explicit dependence on the shower energy. 

The impact of the corrections is clearly illustrated in Figure [HI left, where the y scan 
of ECAL is shown for the raw and, respectively the corrected data. The low energy tail of 
the energy distribution is greatly reduced (Figure [HI right). The resolution loss when going 
from the out of gap events to all the events (with corrections applied on the energy) is of 
the order of 10%. The corrections do not degrade the linearity. 

When tracking information is available, it is possible to precisely calculate the shower 
position within each wafer. Subsequently, the ratio of the active to non- active areas crossed 



before (black triangles) and after cor- 
total energy as function of the shower 



barycentre y (left)^ and energy distribution (right). 
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by the shower according to a mean shower shape can be estimated and used to correct the 
energy recorded in each layer. 

7 Shower development 
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Figure 9: Longitudinal shower profile for the 
data (points with statistical uncertainties) 
and Monte Carlo simulation (histogram). 
The smooth curve is the used parametrisa- 
tion of the shower profile. 



Figure 10: Left: distribution of radii for 90% 
and 95% signal containment (3 GeV e~, nor- 
mal incidence, centre of wafers) . Right: radii 
for 90% and 95% signal containment of e~ 
as function of the beam energy. 



Only events outside the interwafer gaps were used to characterise the longitudinal de- 
velopment of the shower. The mean energy distribution is well fitted by the standard 
parametrisation, ^{t) = c exp(— ^t), where t is the calorimeter depth, c is an overall nor- 
malisation, a and /3 are constants (Figure [9]). The position of the shower maximum grows 
logarithmically with the beam energy. 

An important issue in the development of a calorimeter is to achieve the smallest possible 
effective Moliere radius, in order to provide the best shower separation. It requires the use 
of an absorber with a small intrinsic Moliere radius {Rm)^ but also the minimisation of the 
gaps between the absorber layers. 

Figure [TOl left shows the event distribution for 90% and 95% levels of signal containment 
with respect to the radius. The results for the various energies studied are summarised in 
Figure [TOl right. The points correspond to the peak position of each radius distribution. At 
90% (95%) shower containment the corresponding radius, often quoted as 1 Rm^ is about 
20 (28) mm. 

The geometry of the ECAL prototype, with 2.2 mm thick interlayer gaps leads to an 
effective Moliere radius which is expected to be larger by a factor of about 2 with respect 
to Rm of solid tungsten. The results from the test beam studies are therefore in agreement 
with expectations. R&D effort towards the use of Si pads with integrated readout is under 
way and will hopefully lead to a significant decrease of the interlayer gap and therefore of 
the ECAL effective Moliere radius. 



8 Spatial and angular resolution of ECAL 

The spatial and angular resolution of the ECAL are studied with the DESY data at normal 
incidence. The shower direction and position at the ECAL front face are constructed on 
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an event-by-event basis using a linear two-parameter chi-square fit to the shower barycentre 
positions in each layer for the x and y coordinates separately. The correlation matrix is 
determined from simulations for each beam energy. 

The fit results are 
compared with the po- 
sition and angle mea- 
sured by the tracking 
system. The expected 
e~ position and direc- 
tion at ECAL front face 
is obtained from a lin- 
ear fit of the drift cham- 
bers. Sources of sys- 
tematic uncertainties as 
residual misalignement, 
material modelling, and 
background rate are esti- 
mated for the extrapola- 
tion to the ECAL front 
face. 

The ECAL resolution, deconvoluted from the tracking errors, is displayed in Figure [TTl 
9 Conclusion 

The Si/W ECAL prototype was presented, as well as the results of the first beam tests at 
DESY and CERN during 2006. The prototype calorimeter was further exposed to beam at 
CERN in summer 2007. Analysis of the data collected is in progress. 
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Figure 11: Resolutions in position (left) and angle (right) as 
a function of the beam energy. The data are shown as points 
with error bars. The simulation expectation is shown by the 
continuous line. 
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